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A empirical model is proposed for the observed variation of the London penetration depth of cuprate 
superconductors with applied magnetic field H below the lower critical field Hc\. It is suggested 
that the Doppler shift of the quasi-particle energies near the d-wave nodes could cause pair-breaking 
and that this could be responsible for the non-linear Meissner effect (NLME) investigated in recent 
experiments. This picture gives a somewhat better description of the observed temperature and 
field dependencies than the original theory proposed by Yip and Sauls. A calculation based on 
weak-coupling BCS theory suggests that when there is a finite supercurrent, the state for which 
pairs are broken in a small angular range around the nodes can have a lower free energy than the 
d^2_y2 state. 



I. INTRODUCTION 

There is now a large body of experimental evidence 
[0-^ showing that the order parameter of all cuprate 
superconductors has unconventional dx2_y2 symmetry. 
The superconducting gap parameter [Q varies as 
A cos 29 in the fcj;, ky plane of reciprocal (fc) space, where 
9 = arctan ky/k^, giving rise to line nodes at 45° to the x 
and y symmetry axes, at which the gap goes to zero and 
the order parameter changes sign. In a seminal paper, 
Yip and Sauls predicted ||^ that there should be a non- 
linear Meissner effect (NLME) arising from the anomalous 
back flow of Bogoliubov quasi-particles near these nodes. 
They showed that at low temperatures (T) the London 
penetration depth A should be increased significantly by 
small applied fields H below the lower critical field Hd, 
and that there should also be non-linear effects in the 
transverse magnetisation 11 . The NLME has been stud- 
ied theoretically by several groups |^-^ who all agree 
that A should increase as at very low temperatures, 
when the Doppler or energy shift d in the quasiparti- 
cle spectrum caused by the flow of supercurrent is larger 
that fc^T. For YBa2Cu306+2; crystals this condition is 
only satisfied below approximately 2 K for fields of or- 
der 0.5Hci or 300 Gauss. At higher temperatures, where 
ksT ^ d, the theory predicts that the non-linear effect 
in A should be much smaller and vary as . In contrast 
to these predictions, several experimental studies 10-1^] 
show that A increases as over most of the temperature 
range studied (typically 1.4 to 80 K), and that dX/dH 
becomes larger at higher temperature. Moreover the ex- 
pected anisotropy in SX{H), an increase by a factor \/2 
when H is applied along the nodal directions, and the 
predicted non-linear transverse magnetic moment have 
not been observed experimentally p^ , p^ . One possible 
explanation for these discrepancies is the occurrence of 
non-local effects 1l7| at low T. However in one of the 



commonly used experimental geometries, where the cur- 
rents and fields are in the Cu02 planes, these effects are 
governed by the c-axis coherence length and would gener- 
ally be small, although it has been argued [|l^ that they 
would still be significant in YBCO. 

It is generally agreed that the experiments are diffi- 
cult and not very reproducible, and it has been sug- 
gested that the lack of reproducibility could be caused by 
trapped flux Q or the sharp edges of the crystals . 
My original motivation was to see whether the unusual 
dependence of A mentioned above and perhaps some 
of the apparent lack of reproducibility could be linked 
with other physical properties which are extremely de- 
pendent on the hole concentration (p). For example the 
normal state pseudogap develops suddenly when p is less 
than 0.19 holes per planar Cu atom | |l^ , [T9t , and this 
could conceivably lead to unusual effects in X{H). The 
thermoelectric power is often large [ pO[ , implying sub- 
stantial electron-hole asymmetry which might also affect 



the NLME. However simple model calculations |21 1 based 
on standard theory seem to indicate that none of these 
mechanisms can cause A to increase linearly with H at 
higher temperatures. 

As discussed in the present paper, a possible expla- 
nation of the experimental findings is that Cooper pairs 
near the nodes of any superconductor with higher order 
angular momentum pairing are actually broken when the 
energy difference between the {k + s) 1 and (— fc + s) I 
states forming a pair with drift momentum 2hs exceeds 
the pair binding energy 2|Aj|. An empirical model along 
these lines gives a reasonably good account of the avail- 
able experimental data, although further experiments on 
the same crystal at closely spaced doping levels may be 
worthwhile. The model gives effects of the same magni- 
tude as Yip and Sauls predicted at low T, but they persist 
to higher temperatures. It could also be checked exper- 
imentally by measuring zinc doped crystals where the 
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NLME is predicted to be shifted to higher fields because 
of increased pair-breaking in zero field [^,^. Exten- 
sions of the present approach to higher Doppler shifts {d 
greater than ksT) could also be relevant for understand- 
ing some properties of cuprate superconductors at higher 
magnetic fields, for example the penetration depth J24] 
and the thermal conductivity in the vortex state pM. 
However its relation to the usual theory involving the 
Doppler shift |^,^ needs a more careful study because 
both pictures predict similar numbers and types of quasi- 
particles at high fields and the Doppler shift approach 
does appear to agree with recent thermal conductivity 
experiments [|28| . 

II. THE EMPIRICAL MODEL 

The total supercurrent density j in a superconductor is 
usually expressed as the sum of the contribution js from 
Cooper pairs with a drift velocity Vs and the contribution 
from the backflow of quasiparticles jqp 0. For a cylin- 
drical Fermi surface(FS), with quadratic in-plane energy 
dispersion, and Vs in the conducting planes, js = nevg 
where n is the total number of carriers per unit volume. 
As far as I know, published theoretical treatments of the 
NLME only consider non-linearity in the quasiparticle 
term, and implicitly assume that the gap parameter 
retains its usual d-wave form. For example 0| , for a cylin- 
drical FS, Equations |^ and || together with the relation 
j = js + jqp, lead to all the non-linear effects proposed 
until now. 

4en de 

/ d9vF cos6l/(v/e2 + A2(6') + mvpVs cos 6) (1) 
Jo 

In this Equation, / is the Fermi function, e the elec- 
tronic charge, vp the Fermi velocity, e and m the energy 
and effective mass of the carriers respectively. The super- 
fluid velocity Vs is directed along the crystallographic x 
axis and the angular dependent superconducting energy 
gap is given by: 

|A,,(r)| = A(T,0) = A(T)|cos20| (2) 

At low T when the Doppler energy shift of one elec- 
tron in a Cooper pair, d — ^m[(vp + VsY — (vf)'^] = 
mvp.Vs ^ ksT, Equations^ and|^ give the non-analytic 
variation A cx \vs\, i.e. A cx \H\ that was predicted by 
Yip and Sauls ||^. In the opposite limit, d ^ fcsT, S\/\ 
is expected to vary as H^/T. Here I suggest that 
Equation ^ does not hold exactly except when Vs — > 0. 
Namely, as sketched in Figure 1(a), in the parts of the FS 
where \d\ > |A^(r)|, A(r) = 0, the Cooper pairs are ac- 
tually broken and do not contribute to the supercurrent 
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FIG. 1. (a) Suggested angular variation of |A(6)| in the 
presence of a finite supercurrent (solid line), (b) Section of 
a cylindrical Fermi surface showing the angular regions of 
width 56 where pairing is suppressed. The loss in current 
density from Cooper pairs in the angular range 56 is equal to 
2evF cos 6 (the current change on moving one carrier from R' 
to R) multiplied by the number of carriers per unit volume in 
the shaded areas R or R'. See text. 

[ pof . This differs from the usual picture, shown in Fig- 
ure 2, where the FS is still displaced near the nodes and 
the pair current there is completely compensated by the 
back-flow when (a) d 3> fcsT and d> A{9), giving strong 
non-linearity (|] in A, or (b) d < fc^T and fc^T > 3A(6') 
and the non-linearity in A is weak Pair-breaking 
is plausible physically because when \d\ > |Aj(T)| the 
energy difference between the two electron states which 
form a Cooper pair is larger than the pair binding en- 
ergy, 2|A^(T)|, and in an s-wave superconductor, where 
A^ is uniform, this condition corresponds to the critical 
(depairing) current. 

To some extent it can be justifled within weak-coupling 
BCS [0 theory as follows. Using the standard d-wave 
pairing interaction |^l|], Vkk' = —Vcos29kCos29k', the 
BCS Equations |Q 3.15 and 3.16 can be used to compute 
the increase in free energy F when A{9,T) is set equal 
to zero in a small angular region 69 near the four nodes. 
For the stationary superfluid this costs an energy: 

SF = (3A{T, TT/2f69^N{0) (3) 

where N(0) is the density of states at the FS for one spin 
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FIG. 2. Partial density of states(DOS) for Doppler shifted 
Bogoliubov quasi-particles at an angle 9. At low T, e.g. 
ksT = 0.05d where d is the Doppler shift, the backflow of 
quasi-particles in angular regions where the superconducting 
gap A{6, T) < d exactly compensates the flow of Cooper pairs 
arising from the shifted DOS in those regions j^. At high T, 
e.g. ksT = d, this does not happen. 

direction, and /3 is a weakly T dependent constant, which 
ranges from 1.49 at T/T^ = 0.16, to 1.19 at T/T^ = 0.67. 
If only two nodes are affected (e.g. for current flow in the 
nodal directions) then (3 is reduced by a factor of two. For 
the moving superfluid, there are no Cooper pairs in this 
small angular region, so the kinetic energy K is reduced, 
and in fact SK varies linearly with S9. However there 
is also a magnetic self energy term (M) arising from the 
magnetic field produced by the supercurrents which has 
the form Q j.A, where A is the magnetic vector poten- 
tial. In the London gauge |0] ^ = mvs/e and I argue be- 
low that jqp is not altered by pair-breaking at the nodes. 
Hence there is a gain in magnetic energy of the form 
mSjs.Vs/e arising from the absence of Cooper pairs near 
the nodes. Referring to Figure 1(b), Sjs-A is given by 
2evF-A, multiplied by the total number of carriers in the 
region R, namely s cosOkpSOn/irkp — nVsSOcosO/nvp. 
Hence for one pair of opposite nodes: 



SM = —nmvl 59 cos^ 6 



(4) 



By integrating over the whole of the FS, from 9 — to 
TT, it can be seen that M = nravl = 2K. Minimising the 
sum of 5F and 5{M -\- K) with respect to 59 gives: 



= 0.64 



mvFVs 



(5) 



when Vg is along a nodal direction. 59 is a factor of 
\/2 smaller when Vg \\ x. This expression has the same 
form and the same anisotropy as the qualitative argu- 
ment given above, the presence of the extra factor 0.64 
worsens the agreement with experiment slightly and has 
been omitted in subsequent discussion. This may be an 



indication that there is another pair-breaking mechanism 
or another term in the free energy varying linearly with 
59^ not covered by the above treatment, or that the num- 
bers are changed slightly by strong coupling effects since 
it is known that for many cuprates A(0) ~ iksTc rather 
than the weak couphng d-wave value, A(0) = 2AAkBTc. 

The effect of the broken pairs on the penetration depth 
can be estimated in a similar way. A crucial question 
is whether the carriers produced by the broken pairs 
will have zero drift velocity in the laboratory frame or 
whether they will be in equilibrium with the pair cur- 
rent and the Bogoliubov quasiparticles on the remain- 
der of the FS. It seems that only in the latter case will 
the present model account for the behaviour of A(H). 
From Figure 2 it can be seen that in both cases (a) when 
ksT < d and (b) when kgT > d and keT > A{T,9), 
the pair supercurrent in that angular region will be bal- 
anced by the backflow, to give almost zero net super- 
current, therefore producing stationary carriers by pair- 
breaking will have almost no effect on A. This is con- 
firmed by numerical calculations based on Equation 
namely, if the contributions to jqp and js from the an- 
gular region 59 are simply omitted there is almost no 
change in A for fc^T ^ d. On the other hand if the car- 
riers from broken pairs have the same drift velocity as 
the Bogoliubov quasiparticles then jqp can be calculated 
from Equation simply by setting A (6*) — near the 
nodes. For ksT ^ d, 5 jqp only varies weakly with d, as 
d^, while 5js varies as d^. Assuming then that the car- 
riers do have the same drift velocity as the Bogoliubov 
quasiparticles, the loss in supercurrent associated with 
the broken pairs is given by 5js = 2nevs59 cos^ 9/Tr for a 
pair of nodes. This is the same formula used above when 
calculating the magnetic self energy of the supercurrents. 
Hence from this formula, 5n/n = 5js/js = —259/n for 
both Vs parallel to the nodal directions (2 nodes with 
cos, 9 = 1) and to the anti-nodal directions (4 nodes with 
cos 61 = 1/V2). 



The penetration depth A(r) = ^mc^/47re^(n — 7 
and neglecting non- linear effects on Uqp = \jqp/evs 
the reasons given above, this changes by: 

5\ _ I 5n 
T ^ 



1-qp), 

I for 



(6) 



2 (n - Uqp) 

Using the relation H — 4ttI/c between the surface field 
(H) and the surface current / per unit length yields H = 
AnjX/c, which together with the above expressions for A 
and 59 (but omitting the factor of 0.64 in Equation |^) 
gives: 

5X{H) _ -f{T)d _evF\H\ \{T)-f{T) 



A(T) 



V2ttA{T) cV27rA{T) 

7t\h\ e(r)A(T)7(r) 



aj(T)- 



H 



(7) 



y2$o '"^'H,{T) 
This expression refers to the case where the current 
flow is at 45° to the nodes. According to the present 
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treatment the effect will be a factor of -\/2 larger when the 
current flows in the nodal directions, as in the work of Yip 
and Sauls $0 is the flux quantum for pairs, ^(T) = 
?iUi?/7rA(T) is a measure of the in-plane coherence length 
for a d-wave superconductor and 7(T) = A(T)^/A(0)^. 
Hc{T) is the thermodynamic critical fleld, if the stan- 
dard Ginzburg-Landau expression for an s-wave super- 
conductor, HciT) = $o/2A/27rA(T)^(r) is used, the con- 
stant a = 0.25. Estimates using parameters for a cylin- 
drical FS give a — 0.16 for the weak coupling d-wavc 
case, (for which A(0) = 2.UkBTc and Hc{0)'^/8tt = 
0.125Af(0)A(0)2) and a= 0.129 for the more realistic case 
A(0) = afcsTc and H H^{0)'^/8tt = 0.089iV(0)A(0)2. 
Equation 7 has similar form and similar magnitude to the 
Yip-Sauls result at low T, but it has a different origin, 
in the present work it arises because the superconduct- 
ing gap parameter has been suppressed near the nodes, 
where it no longer obeys Equation ^j. The extra factor 
A(r)^/A(0)^ results from the assumption that jqp, i.e. 
Uqp, is not changed at all by the presence of pair-breaking. 
If it were changed in the same ratio as n, i.e. if Snqp/uqp 
= Sn/n, then the A(r)^/A(0)^ factor would be absent and 
the temperature dependencies would be in much better 
agreement with experiment. But there does not seem to 
be any theoretical justification for such an "ad hoc" hy- 
pothesis. The behaviour in Equation 7 will be cut 
off at low fields by the presence of pair-breaking from 
other sources, for example at low T from impurities in 
the Cu02 planes and at higher T by inelastic scattering 
processes. This effect can be taken into account by re- 
placing |iy I by the empirical expression y/H^ + H^ — Hq, 
where Hq is a measure of the pair-breaking effects in the 
absence of an applied field. 



III. COMPARISON WITH EXPERIMENT 

Experimental data |jl^ for the field dependence of Aq 
for a detwinned YBa2Cu3 06.95 crystal are shown in Fig- 
ure 3a, together with fits to Equation 7, including the 
cut-off parameter Hq. The linear terms in these fits cor- 
respond to d\a/dH values of 0.033 ± 0.001 A°/Gauss at 
both 4.2 and 7 K. In contrast to the Yip-Sauls result, 
the slope calculated from Equation 6 is T-independent 
at these temperatures. Using the experimental value fl^ ] 
Aa(T = 0) = 1600 A, leaves one free parameter ^a(O) 
which has to be equal to 16.4 Ato account for the data in 
Figure 2(a). This value is quite compatible with £,a^b{0) 
= 11 Aobtained previously by analysing the specific 
heat above Tc in terms of Gaussian fluctuations bearing 
in mind that ^a-b{0) contains a contribution from the Cu- 
O chains, depends on the FS average of vp and is related 
to the RMS average of A(0). Taking 2A(0) = eksTc, the 
value of ^a(O) obtained in the present work corresponds 
to an acceptable value for vp = 1.9 x 10^ cm/s in the 
nodal directions. Using Equation 6 and the same values 
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FIG. 3. (a) Field dependence of the a-axis penetration 
depth of YBa2Cu3 06.95 measured for a detwinned crys- 
tal at 4.2 K (circles) and 7 K (squares). The dashed and 
solid lines are fits to the data using Equation 7 together with 
cut-off fields Hq of 23 and 42 Gauss respectively, (b) Mea- 
sured temperature dependence of dAab/d// (squares) for 
a YBa2Cu3 06.95 crystal with Tc = 91.4 K. The dashed lines 
show the behaviour calculated according to Equation 7 for 
various magnetic fields corresponding to d/\/2A(-|,T = 0) = 
0.001 to 0.004. The solid line shows the effect of omitting the 
A(T)^/A(0)^ factor, see text. 

of A(0) and Aa(0), the value of Hq = 2,0 ± 10 Gauss, 
obtained from the fits, gives a very small pair-breaking 
energy of 0.5 ± 0.16 K. This compares well with typical 
values of the pair-breaking temperature T* < 1 K ob- 
tained from the T dependence of A in zero field [ p"3p^ ]. 

The measured T dependence Q of dAab / dH is shown 
in Figure 3b. The fit to Equation 7 using the d-wave 
behaviour of A(r) and A(T), calculated in the weak- 
coupling limit for a cylindrical FS, gives the correct gen- 
eral behaviour with some deviations at high T. Exact 
agreement is not to be expected in view of the fact that 
the weak-coupling ratio 2A(0) — 4:.28kBTc is smaller than 
typical experimental values, 2A(0) = QksTc- However as 
also shown in Figure 3b, omission of the A(r)^/A(0)^ fac- 
tor gives much better agreement with experiment. 

Experimental data for SXab{H,T)/ Xab{0,T) versus 
H/Hc{T) from the work of Maeda and colleagues [fo 11 
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FIG. 4. (a) and (b) Experimental data for the in-plane pen- 
etration depth of a Bi:2212 crystal as a function of the 
applied field, H || c, normalised to the thermodynamic critical 
field Hc{T) at various temperatures. The solid lines show fits 
to Equation 7. (c) Values of the parameters a (left-hand scale) 
and Hq/Hc{T) (right- hand scale) determined from these fits, 
omitting the A(r)^/A(0)^ factor, which as shown by the solid 
line gives the wrong T dependence. The dashed line shows 
the T dependence expected if the pair-breaking parameter Ho 
is constant. 



for Bi:2212 crystals are shown in Figure 4(a). Hc{T) was 
derived experimentally pO| , pT[ | for the same crystals by as- 
suming that the lowest field H* , at which some hysteresis 
could be detected in field sweeps, was a measure of the 
lower critical field Hd {T) . Then He [T) was obtained by 
taking a plausible value k = 100 for the Ginzburg-Landau 
parameter. It can be seen that the data can be fitted to 
Equation 7 at all reduced temperatures [t — T/Tc) stud- 
ied. However as shown in Figure 4(b) the values of the 
slopes are essentially constant rather than increasing as 
A(r)^/A(0)^ as predicted by Equation 7. Again this ex- 
tra factor of A(T)^/A(0)^ is a weak point of the present 
model. Even without the latter factor the magnitude 
of the slopes, i.e., a, are about a factor of 10 larger than 
that calculated above. It was suggested Q that this was 
mainly because the effect of the geometrical barrier HJ] 
leads to Hci{T) and hence Hc{T) being overestimated. 
The crystals used were slightly overdoped with Tc val- 
ues of 90 K. Although it is known that the condensation 
energy of Bi:2212 can be reduced by a factor of 2-5 by 
the presence of the normal state pseudogap Q, even 



for such high values of Tc, this effect is probably taken 
into account by using experimentally determined values 
of He{T). The pair-breaking parameter Ho{T)/Hc{T) 
obtained from these data is also plotted in Figure 4 (b). 
Hq is essentially constant in the temperature range stud- 
ied, showing that inelastic scattering processes are not 
strong enough to cause pair-breaking. 



IV. SUMMARY AND CONCLUSIONS 

I have argued that the non-linear Meissner effect in 
the cuprates can be understood in terms of an empiri- 
cal model in which the flow of a finite supercurrent de- 
stroys superconducting order in a small angular range 
near the c?-wave nodes. This model accounts for some of 
the experimental data and would mask the unusual tem- 
perature dependence predicted by Yip and Sauls. The 
general viewpoint outlined here is reminiscent of pub- 
lished theoretical work suggesting that application of 
high fields could actually alter the symmetry of the or- 
der parameter, for example by introducing a d + id [3^ ] 
or a d -|- idxy component | |36[| . However for experiments 
at low fields it seems that the empirical pair-breaking 
picture described here could be adequate. The present 
approach, which suggests that the pair-breaking is an in- 
trinsic property, still predicts anisotropic effects for which 
there is presently no experimental support. Although 
in this regard one should note that there are potential 
complications in separating CU-O2 plane and Cu-0 chain 
contributions to the penetration depth in YBCO and that 
in the present model the details of any non-linear effects 
in the magnetisation will be altered by the different be- 
haviour of the free energy. The present work supports 
the empirical conclusion ||lo|-|l^ that even though the 
high temperature NLME is not predicted by the origi- 
nal Yip-Sauls theory, it is nevertheless a characteristic 
feature of 'non s-wave' superconductivity. 

There are indications that measurements of the non- 
linear Meissner effect as a function of hole concentration 
could be used as a means of investigating the effect of 
the normal state gap on the superconducting properties 
and for investigating pair-breaking mechanisms in zero 
applied field. 
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